#if 1
#define MECAFLUIDE_CL(TIME)                                             \
    /*addCL<cl::paroi_mobile,cl::ParoiB>(vf::cst(0.))                   \
      + */                                                              \
    addCL<cl::paroi_mobile,cl::Paroi>(vf::vec(vf::cst(0.),vf::cst(0.) ) ) \
    /*    +                                                             \
          addCL<cl::dirichlet_vec,cl::Paroi>(vf::vec(vf::cst(0.),vf::cst(0.) ) )*/ \
    +                                                                   \
    addCL<cl::neumann_scal,cl::Entree>(-vf::cst(1.e4)*( 1-vf::cos(M_PI*TIME/(2.5e-3)) )*( vf::chi(TIME< 0.0025)) ) \
    /*addCL<cl::neumann_scal,cl::Entree>(-vf::cst(1.332*1.e4)*( vf::chi(TIME< 0.001)) )*/ \
    +                                                                   \
    addCL<cl::neumann_scal,cl::Sortie>(vf::cst(0.))                     \
    /**/
#else
#define MECAFLUIDE_CL(TIME)                                             \
    addCL<cl::paroi_mobile,cl::ParoiB>(vf::vec(vf::cst(0.),vf::cst(0.))) \
    +                                                                   \
    addCL<cl::paroi_mobile,cl::ParoiH>(vf::vec(vf::cst(0.),vf::cst(0.))) \
    +                                                                   \
    addCL<cl::dirichlet_vec,cl::Entree>( 343.99*(0.5-vf::Py())*(0.5+vf::Py())*(- 1357*vf::pow(TIME,9 ) \
                                                                       + 7443*vf::pow(TIME,8 ) \
                                                                       - 17099*vf::pow(TIME,7 ) \
                                                                       + 21255*vf::pow(TIME,6 ) \
                                                                       - 15356*vf::pow(TIME,5 ) \
                                                                       + 6379*vf::pow(TIME,4 ) \
                                                                       - 1368*vf::pow(TIME,3 ) \
                                                                       + 97*vf::pow(TIME,2 ) \
                                                                       + 6*TIME)*vf::oneX() ) \
    +                                                                   \
    addCL<cl::neumann_scal,cl::Sortie>(vf::cst(0.))                     \
    /**/
#endif


#define MECAFLUIDE_USE_LAGMULT 0                \
    /**/

#define MECAFLUIDE_FORCE_VOLUMIQUE(NU)                                  \
    vf::vec(vf::cst(0.),vf::cst(0.) )                                   \
    /**/
